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Abstract 

In high energy heavy ion colhsions the pion multipUcity is large, and one 
might expect that pions are radiated semi-classically. The axial symmetry of 
the collision and approximately zero isotopic spin of the colliding nuclei result 
then in peculiar isotopic spin - azimuth correlations of the produced pions. 
These correlations are easy to test - and should be tested 
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Central collisions of heavy ions at high energies are expected to produce 
many pions (see for instance so many that the typical number of pion 
quanta per unit cell of phase space may substantially exceed unity. Under such 
circumstances, one may speculate about the possibility that, in the early stages 
of the collisions, a classical pion field develops, before subsequently decaying 
into pions (2|, ^, ^. One interesting possibility showing how this could happen 
following chiral symmetry restoration has been investigated recently |^]. Inde- 
pendently of the detailed dynamics, the formation of an intermediate classical 
pion field could result in certain correlations in the production of pions with 
definite isospin components. The main point of this letter is to show that, 
indeed, rather peculiar correlations can be indicated from symmetry consider- 
ations only - the role of dynamics is to guarantee that the pion radiation is 
so strong that it can be treated semi-classically; testing the correlations is a 
way to check that hypothesis. These correlations, quite distinct from the much 
studied Bose Einstein correlations, find their origin in the special correlations 
which develop between spatial and isospin coordinates in the solution of the 
non linear field equations obeyed by classical pion fields. 

Mathematically, one can write the amplitude of N pions' production with 
the help of the Lehmann-Symanzik-Zimmermann formula, where the A^-point 
pion Green function is presented through the functional integral over the pion 
fields: 

A^'-^^ih ...%) = ^^lim ^ J Dtt" J DJ" W[J] exp (^iS[7r] + i J d'^xn^J" 

N 

■ n / d'xne^'-^-i-dl - ml)7r^-{xr.). (1) 

Here ^[vr] is the effective pion action and J is the source formed by the col- 
liding nuclei; one has to integrate over the sources with some weight functional 
W[J] which summarizes the dynamics of the collision. If the source J is in a 
sense "large" , the functional integral can be evaluated in the saddle-point ap- 
proximation, where the saddle-point pion field is the solution of the equations 
of motion, 

+r{x) = o, (2) 



supplemented with the radiation condition at large distances, which at — > 
reads 
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Let us denote the solution of these eqs. as 7r'^(r, t). In the leading WKB 
approximation one replaces the pion fields everywhere in eq. dll) by this saddle- 
point field. In the next approximation quantum fluctuations about the classical 
field 7r"(r, should be also taken into account. We do not discuss dynamical 
questions in this letter but concentrate solely on general symmetry considera- 
tions. We use two assumptions: 

1) the collision is axially-symmetrical; 

2) the collision is isotopically-symmetrical; that would be an exact state- 
ment in case of the zero isospin of the colliding nuclei and is anyhow a good 
approximation for pions produced in the central rapidity region since a transfer 
of any quantum numbers over a long rapidity range Ay is suppressed at least 
as {Ay)-\ 

The first assumption means that the saddle-point pion radiation field can 
be sought in the axially-symmetrical "hedgehog" form: 

7r'^(r,t) = O'^* {n^,P^{p,z,t)+n\\iP\\{p,z,t)) (4) 

where P±^\\ are functions of the distance p from the beam axis, of the dis- 
tance z from the collison point and of time t; n^\(±) are unit vectors parallel 
(perpendicular) to the beam axis and is an arbitrary 3x3 orthogonal ma- 
trix. (In principle, axial symmetry does not contradict higher harmonics in 
the transverse plane with n_L replaced by a unit vector with the components 
{cos m(f), sinm(/>) where (j) is the azimuthal angle and m is an integer. We 
shall work with m = 1, i.e. assume that n_L points in the radial direction, and 
introduce the arbitrariness in the choice of m only in our final results). The 
type of correlations between space and isospin degrees of freedom is a familiar 
feature of solutions of non linear field equations satisfied by the pion fields, as 
illustrated for example in the Skyrme model (see, e.g. [Q). The second assump- 
tion means that the saddle point is degenerate in the global isospin rotation 
0°'^, so that the functional integral in eq. (||) reduces to the integration over 
all possible orientations of the pion field in the isospin space. The action and 
the source-weight functionals are taken at the saddle-point values of J*^ and tt"" 
and provide an overall normalization factor ^/J\^ which may be a function of 
the total 4-momentum of the pions P^, but now we are not interested in this 
factor. 
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At large distances / time the isospin source J*^ dies out and one can also 
neglect the non-linearity of the pion effective action. Therefore, eq. (Q) reduces 
to the free Klein-Gordon eq. at large distances, and we are guaranteed that 
the Fourier transform of eq. (^) has a pole at k^ — k'^ — = ni^, which cancels 
out in the LSZ leg amputation procedure (see eq. ([l|)). We have therefore: 

lim / d^xe'''''{-dl - ml)TT''{x) = O"' (^±iF±{k±,ku) +'kuiF»{k±,h\ 

= O'^'Fiik) (5) 

where -Fj_j| {k±, k^^ ) are related through Fourier transformation to P±^\\ of eq. (Q). 
Since P±,\\ are real the functions -Fi(k) are purely imaginary, with -F*(k) = 
Fi{-k) = -Fi{k). 

Squaring the amplitudes, summing over the isospin a = 1, 2, 3 of the pions 
and multiplying by the phase space factor, one gets for the N pion production 
cross-section: 

"^""^ = ^ /'^^^'^^^ / n ^,2.6^{kl-ml){2.r 

n=l ^ ' 

• OfFi(k„)02"'F;(kn)<5(^)(^/. - SV), (6) 

where knfj, are individual 4-momenta of the produced pions, Oi^2 are isospin 
orientation of the pions in the amplitude and the conjugate amplitude, respec- 
tively, P is the total 4-momentum of the produced pions; the factorial accounts 
for the identical particles and 6+{k'^ — m^) = 6{k'^ — m^)6{ko). 
Writing the 4-momentum conservation restriction as 

(27r)^5W(P -Y,kn) = I fi4i?e*(^-^)-^i:('="-^), (7) 



we get factorized integrals over the momenta of produced pions: 

J ^^2n6+{e -ml)e-^(>^-^^F,{k)F*{k) =T.,{Ro,R). (8) 

This tensor is further on contracted with the relative orientation matrix 
0^2 = 0"*02"' = (0^02)*-'. We note that in integrating over the SO {3) rota- 
tions one can use the Haar measure property, dO = d{CO) = d{OC), where C 
is an arbitrary orthogonal matrix, so that 

J J dOid02 = J J dOidOi2 ; J dO = I. (9) 
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The total cross section being a sum of a^^^ over becomes thus a series 
for an exponent, and we get 

^tot =M J dOi2 J d^Rew [i{P • R) + Tij{Ro,R.)0'^2\ ■ (10) 

Let us use the frame in which the total momentum of the produced pions 
is zero, P = 0, Pq = E where E is the total energy of the pions. Since both E 
and the function J^, proportional to the probability of the pion production, are 
presumably large, one can integrate over Rq, R and Ou by the saddle-point 
method. Let us parametrize the relative orientation matrix O12 in terms of a 
unit 4-vector = 1, = (1, —ir)^ , 

O12 = i^Tr{u^T^T''uyTlT^) = (1 - 2v?)5ij + 2uiUj + 2uQUkeijk- (H) 

We expect the saddle point to be at 

ussO, RssO, i?o ~ 4- (12) 

E 

Expanding the exponent in eq. (|lO|) around this point, and using the explicit 
form of J- given by eq. (|8|), it can be shown that, at small Rq, the saddle- 
point condition is indeed satisfied. Without knowing the explicit form of the 
functions -Fi(k) we cannot prove that there are no other saddle-points, but we 
shall disregard such possibilities. Note that the maximum at u = corresponds 
to 0^/2 = "^'^ i-e- to the case when the pion isospin orientation in the conjugate 
amplitude is the same as in the direct amplitude - not an unnatural result. 

We thus write the total cross section as 

cj*°* j dRQ exp [i{ERQ) + TuiR^, 0)] . (13) 

In what follows we shall use this relation to remove the unknown normalization 
factor J\f. 

We next turn to the 1,2,... particle inclusive cross sections. All of them are 
directly derived from eq. (P) where one skips integration over 1,2,... momenta 
and summation over the 1,2,... isotopic subscripts. Thus, the 1-particle inclu- 
sive cross section is given by (we use the abbreviation (dk) = (i^k/2A;o(27r)^): 

^ = ui^^^""^ J dO,dO,Of Or^F.(k)i^nk) 
•^E Tlv^W / I[idK)iFiK)OuF^iK))i2n)'6^'\P-k-Y^K). (14) 

A''=l ^ ' n 
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(A) 

Here u, are the isospin "polarization" vectors, 



/0\ 


V 1 / 



u 



1 

7^ 



1 

7^ 



/ 1 \ 



V J 



(15) 



b ^ ^ b ^ ^ b 

to be contracted with the isospin orientation matrices Oi^2- The superscript A 
refers to the isospin component of the observed pion; there is no summation on 
A. 

The 2-particle inclusive cross section for production of the pion of sort Ai 
with momentum ki and of the pion of sort A2 with momentum k2 is 



. 00 

• J dOidOsOf 02'Fi(ki)F;(ki)Of OiVfc(k2)Fr(k2)AA J] 



I n (dkn) (F(k„)Oi2F*(k„)) (27r)4(5(4) (P-ki-k2-J2 ^n), (16) 



. N-2 



and so on. Writing the 4-momentum conservation 5 function with the help of 
an auxiliary integral as in eq. (|^, we again obtain the exponential series, so 
that 

"^^^^^ ^dOidOi20fof'=o^^F,(k)F;(k) 



idk) 



(A) 



N / (i^iJexp i{P-R) 



R)+Ti,{R)Og 



where we used = Oi^O\^. If the momentum of the observed pion is negli- 
gible as compared to the total momentum P of the pions and if the integration 
over O12 and R is performed about the presumably steep saddle point given by 
eq. ([T2I), we obtain: 



dcr^^'' {A) *(A) 
(dk) ^""^ 



.iA)^*(A)l^bb' , 
3 



H 'V 'T^r 5*^Fi(k)F;(k)a 



dOiOfof^Fi(k)F;(k)cj*°* 
1 



tot 



-F,(k)F/'(k)a*°'. 



(18) 



We see that the inclusive cross section is directly related to the square of the 
Fourier transform of the classical pion radiation field, - a most natural result. 
Also naturally, we find identical cross sections for vr''", vr" and vr*^ production. 
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The average multiplicity is obtained by integrating the inclusive cross section 
over k and summing over A = 7r~^ , 7r~ , vr*^ : 



(N) « J {dk)Fi{k)Ft{k) « Tu{R^, = 0). 



(19) 



It should be kept in mind that if this integral is not convergent by itself at large 
k it should be cut at least at P in accordance with a more precise eq. ([T7|). 
In eq. ( |T7|) we used the following formula for averaging over isospin rotations: 



^b'k 

3 



(20) 



To calculate the 2-particle inclusive cross section we need a formula for averag- 
ing over four matrices: 



1 

30' 
1 

30' 



+ ^S'^5ik^_^bb' ^cc' _ gbcgb'c' ^ 4jfec'^6'c)_ 



(21) 



One can check this formula by applying various contractions and reducing it to 
eq. (|20|); an alternative method is to note that O''* is a Wigner D function for 
isospin 1, and using the Clebsch-Gordan machinery. 

Starting from eq. (^) and repeating the same steps as above we get for the 
2-particle inclusive cross section: 



(dki)(dk2 



(A^) MA^) XA2) AA2) 

H %' "c "c' 



30 



26''^' 6^"' {2V -W) + {6'"'6'''''' + 6"''' 6'''''){-V + 31^) 



be' s;b'c\ 



a 



tot 



tot 



3V + W 
W-2W 
^ 2V + AW 

where we have introduced the abbreviation: 



a 
"30 




(22) 



y = |F,(ki)|2|F,-(k2) 



W = \Fiiki)Fi{k2) 



(23) 



The upper line in eq. (^) (case 1) refers to the vr+vr^, vr^vr^, vr+vr^" or vr^vr^ pro- 
duction, the second line (case 2) refers to four other combinations, vr'^vr'^, vr^Tr*^, 
Tr^7r~^ ,Tr^ it" , while the last line (case 3) corresponds to the production. 
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It should be noted that, if the two observed pions are not identical, they are 
distinguished by the momenta ki,k2. The three possibilities in eq. ( p^ reflect 
three possible isospin states (T = 0,1,2) which can be formed by a pair of 
pions. Since, however, the three cross sections are expressed through only two 
functions, we get a relation: 



2cr 



(1) 



(24) 



Let us now investigate eq. (^21). If one sums up all 9 possible combinations 
of pion pairs, one gets 

do-"" 



|F,(ki)|V:,-(k2)|V*°* 



(25) 



which is independent of the angle between the two pions. Further on, inte- 
grating eq. (p5|) over the momenta ki^2 and recalling eq. ([l^) for the average 
multiplicity, one finds 



{N{N-l)) 



1 

lot 



(dki)(dk2) 



{Ny, 



(26) 



(dki)(dk2) 

which is the dispersion law of the Poisson distribution. It should be stressed 
though that for charged or neutral pions separately there is no Poisson distru- 
bution! Also, we would expect a deviation from the Poisson distribution at 
the end point of the spectrum where the 4-momentum conservation law from a 
more accurate eq. (|l6|) imposes additional correlations. 

The structure denoted as W depends on the azimuthal angle between the 
two pions, 

2 



W 



cos((/)i - (p2)ktktF±{ki)F±{k2) + klkMki)F»{k2 



while 



V 



{ktfFlik,)) + {klrF■^{k^)\ [{ktrFl{k2)) + {kl^YF-^ih 



\2 Tp2 



2 Tp2 , 



^2 772/ 



(27) 



(28) 



Out of the three cross sections mentioned in eq. ( p^ ) one can construct angle- 
dependent and -independent combinations: 

^tot 



2a«+a 



(2) 



2a(2)+a 



(3) 



tot 



-^(ki,k2). 



(29) 
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Eqs.(p^29) summarize our result for double-inclusive pion production. How- 
ever, it might be useful to make a prediction which is independent of the dy- 
namics hidden in the Fourier-transformed pion fields To this end we 
restrict ourselves to pions with zero rapidity, i.e. to the case A^l — — 0? 
that W/V = cos^((/)i — (j)2)- This quantity is obtained by taking the ratio of 
differential double-inclusive cross sections, say, 

4^(1) _ 3^(2) . ^ 

2^(2) +^(3) =^"^ i<Pi-^2). (30) 
Another way to isolate the azimuthal angle dependence is to normalize to the 
single-inclusive cross sections. For example, we predict at /of 2 = 0: 

At this point an experimentalist may derive correlations for his own favourite 
(charged or neutral) pairs of pions. Let us recall finally that the axial symmetry 
does not contradict higher harmonics in the transverse plane, with the replace- 
ment cos(0i — (^2) cosm{(j)i — (f)2) where m is an integer. 

D.D. would like to thank Maxim Polyakov and Michal Praszalowicz for a 
helpful discussion, the Institute for Theoretical Physics-II of the Ruhr Univer- 
sity at Bochum for hospitality during the completion of this paper and the 
A.V.Humboldt - Stiftung for a support. 



References 

[1] H.Satz, Nucl. Phys. A544 (1992) 371c 

[2] A.A. Anselm, Phys. Lett. 217B (1989) 169; A.A. Anselm and M. Ryskin, 
Phys. Lett. 266B (1991) 482 

[3] J.D. Bjorken, Acta Physica Polonica B23 (1992) 561 

[4] J.R Blaizot and A. Krzywicki, Phys. Rev. D46 (1992) 246 

[5] K. Rajagopal and F. Wilczek, "Emergence of Coherent Long Wavelength 
Oscillations After a Quench: Application to QCD" PUPT-1389, lASSNS- 
HEP-93/16 , March 93 

[6] A. P. Balachandran, G. Marmo, B.S. Skagerstam and A. Stern, "Classical 
Topology and Quantum States", World Scientific (1991) 



9 



